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Abstract 

We use a nonpertubative approach which combines coupled channel Lipp- 
mann Schwinger equations with meson-meson potentials provided by the lowest 
order chiral Lagrangian. By means of one parameter, a cut off in the momentum 
of the loop integrals, which results of the order of 1 GeV, we obtain singular- 
ities in the S-wave amplitudes corresponding to the a, fo and ao resonances. 
The 7T7T — ► 7T7T , 7T7T — > KK phase shifts and inelasticities in the T = scalar 
channel are well reproduced as well as the Tr°r] and KK mass distributions in 
the T = 1 channel. Furthermore, the total and partial decay widths of the fo 
and ao resonances are properly reproduced including also the decay into the 77 
channel. The results seem to indicate that chiral symmetry constraints at low 
energy and unitarity in coupled channels is the basic information contained in 
the meson-meson interaction below \fs ~ 1.2 GeV . 
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1 Introduction 



The understanding of the meson-meson interaction in the scalar sector is still prob- 
lematic. There is some debate about the spectrum of hadronic states and even more 
about its nature. Below y/s = 1.2 GeV, which we will study here, the existence of 
a broad scalar-isoscalar meson around 500 MeV has had permanent ups and downs. 
However, the / (980) (also called S*, I G (J PC ) = 0+(0++)), and the a (980) (also 
called 5, I G (J PC ) = 1~(0 ++ )) mesonic states are well established experimentally, al- 
though there is still debate around their decay widths and particularly about their 
nature. Following the discovery of the f and a Q 0, several proposals were made 
about the nature of these states: qq states [ 3-10], multiquark states [JTlj, |12| or KK 



molecules [0, [14], [15], [RJ . Other works argue against the qq nature of the states [17 
It is also interesting to investigate the structure of these states in order to better 
isolate possible candidates for glueballs and other states rich in gluons. For instance, 
in ref. |18|] a glueball state at 992 MeV (5*1(991)) is predicted in addition to the 
/o(980) and the a which have another nature. In ref. |19[] a glueball at energies 
below 0.7 GeV is also predicted which could be a candidate for the o. However, its 
narrow width around 60 MeV is considerably smaller than the 400 MeV associated 
to the conventional a width in the nir interaction. On the other hand in more recent 



calculations p(J the states with a significant component of gluons are only predicted 
at energies around 1.5 GeV or higher. 

The situation in the scalar sector contrasts with the vector and tensor sector 
where the constituent quark models are rather successful in the interpretation of the 
spectrum and properties of the particles. 

Coming back to the meson-meson interaction in the scalar sector, a large fraction 
of the work done consists in a parametrization of the amplitudes respecting general 
principles H §, || |7|, ^TJ while some other works use models inspired in QCD [11], [12 
|13|, |T4"|, |T5[| or phenomenological ones based on the exchange of mesons IT 



Along the lines of the amplitude parametrization it is interesting to quote the work 
of [^, where it is shown that one pole close to a strong threshold (f.i. KK) in 



the II Riemann sheet (we follow the notations of ref. p4|) indicates the presence of 
a molecular state resulting from the forces between the participant mesons (KK, ittt 
in T = or KK, tjti in T = 1). Conversely, the presence of two poles close to 
threshold, one in the II sheet and another one in the IV sheet, would indicate a qq 
state. The analysis in ref. || |6|, 0] favours this latter interpretation, but in ref. 



the first interpretation is advocated, since the pole found in the IV sheet is far from 
threshold. Following this debate, in ref. |7j it is argued that the J/ip — > <p7T7r,<pKK 
decay is crucial in order to distinguish between the two interpretations quoted above, 
questioning the conclusions of [2I| and offering a good description of these data. 
However, in ref. |l(| the fo, appears as a pole in the II sheet and a description of 
the J /if) — > 07T7T, (pKK data of the same quality as in |7[] is obtained. The main 
sources of discrepancy between J7j and [^] come from the different ways in which the 
background is treated. 

Among of the QCD inspired models, in refs. JO], the q 2 q 2 states in a bag 
model is discussed and a rich spectrum of SU (3) multiplets is obtained which would 
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accommodate the a, fo, ao mesons in a same nonet. On the other hand the work in 
ref. |T^, |15| starts from a colour confining Hamiltonian with hyperfine interaction 
and studies the same q 2 q 2 system. According to this latter work the rich spectrum 
predicted in ref. [[II]] disappears because the states separate into colour singlets, 
exception made of two weakly bound KK states with T = and T = 1 which are 
identified as the fo and ao respectively (the authors also warn about the influence 
in such states of other meson-meson components due to the coupling of different 
channels) . 

In ref. [§, |9|] the authors use a unitarized version of the quark model and conclude 
the existence of a a and also that the fo, a are manifestations of ss and qq states 
respectively, although they spend most of their time as KK components. 

In ref. [[16|] the Jiilich meson exchange model [25f is extended to account for the 
meson-meson interaction in the T = and T = 1 channels. A coupled channel 
approach with the KK and nir channel in T = and KK and nt] channel in T = 1 
was followed and the fo and ao states appeared as poles of the t-matrix. 

The different approaches which are successful in reproducing the meson-meson 
scattering amplitudes rely upon a relatively high number of parameters which are 
adjusted to the data, and vary between around 25 in 0, [| ||, 0, ^l|] and 5-6 in |8|, [15 . 

One of the properties which has been used to discriminate among models is the 
77 decay width of the fo, ao states p6fl , and typically has been presented as a support 
for the KK molecule nature of the fo and ao states |27] . However, the experimental 



errors are still relatively high and the data can be accommodated in several models 
A brief summary of the present situation in the scalar sector can be found in ref. [p8 



Another different approach to the meson-meson scattering is chiral perturbation 
theory ||, g(| ||. 

Calculations have been carried out to one loop and one needs counterterms which 
require the use of 10 parameters fitted to experimental data. The approach is valid at 
energies below 600 — 700 MeV but such as the perturbative calculations are carried 
out they show obvious limitations to face the singularities in the t matrix. 

From the former discussion we can see that a key ingredient in most of the models 
which successfully describe the amplitudes in the L — 0, T = 0, 1 sector is the solution 

1331. On 



of coupled channel scattering equations starting from some potential |15|, 
the other hand the success of chiral perturbation theory as an alternative approach 
to the dynamics of QCD, using mesonic instead of quark degrees of freedom, makes 
the use of the chiral Lagrangian extremely appealing. The combination of these two 
factors is the main contribution of the present work. In this sense similar steps in the 
baryon-meson sector have been done in refs. 



34 



| with a remarkable success. We 
shall see that this is also the case here. Our model accounts automatically for unitary 
and analyticity and leads to poles in the scattering matrix for the a, fo and ao states 
below a/s = 1.2 GeV. It predicts the mass and partial decay widths of the fo and 
ao resonances, as well as the different scattering amplitudes, in good agreement with 
experiment and requires the use of only one parameter which is a cut off in the loop 
integrations, following also the spirit of the approach of refs. |34|, |35| . This cut off is 
around 1.2 GeV, very close to the value of the scale of chiral symmetry breaking, A x 
36fl , which plays the similar role as a scale of energies as our cut off. 
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2 L = 0, T = 0, 1 strong amplitudes in lowest order 
of chiral perturbation theory 



We start from the standard chiral Lagrangian in lowest order of chiral pertur- 
bation theory (xPT) C 2 which contains the most general low energy interactions of 
the pseudoscalar meson octet [29-32] in this order. This Lagrangian will provide 
the potentials which will be used in the coupled channel scattering equations. The 
interaction chiral Lagrangian is given by 

where the symbol <> indicates the trace in the flavour space of the 577(3) matrices 
appearing in $ and M, f is the pion decay constant and the matrices $ and M are 
given by 
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where in M we have taken the isospin limit (m u = m d ). From eqs. (1) and (2) we 
can write the tree level amplitudes for KK, irix and nr] (we take r] = rjg) 



a) K + (k)K~(p) -> K + (k')K-(p') 
t a = -^p{s + t-2u + 2m 2 K ) 

b) K°(k)K°(p) -> K°(k')K°(p') 
h = t a 

c) K + (k)K~(p) -> K°(k')K°(p') 
t = H 

d) K + {k)K~{p) -> 7r + (A;')7r-(p / ) 

*d = -q7^(s + t - 2m + m|- + m*) 

e) K + {k)K-{p) -> 7r (fc')7rV) 

* e = -1272 (2s - i - u + 2m^ + 2m;;) 

f) K°(k)K°( P ) ^7T + (k')7T-( P ') 

tf = -qtK-s + m - 2t + + m£) 

g) K°(k)K°(p) -> Tr^W) 

^3 



(3) 
(4) 
(5) 
(6) 
(7) 
(8) 
(9) 
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h) K+(k)K-(p)^n°(k') V (p>) 

t h = -$(2s -t-u+lml- \m\) ^ ] 



i) K°(k)K°(p) -> 7T (k')r](p') 
U = —th 

j) n°(k)r)(p) -> tt°(A;' 

, m| 

r i — — 3/2 



k) 7r°(fc)7r°(p) -> n° (k')ir (p') 

P 



f, - -HH. 



1) 7r+(ifc)7r-(p) ^ 7T°(A;')7rV) 
tj = -3^(2s-w-t + m2) 



(11) 

(12) 
(13) 
(14) 
(15) 



m) 7r + (A;)7r (p) — > 7r + (A; / )7r (p') 
*m = -372(5 + t-2w + 2m 2 r ) 

where s = (A; + p) 2 , t = (k — k') 2 , u = (k — p') 2 

In order to obtain the S-wave amplitudes in the different isospin channels we 
construct the isospin eigenstates projecting over and S-wave satate . We have 

T = 

\KK >= -j=^f(q)\K^q)K-(-q) + K (q)K°(-q) > (16) 
I™ >= -j=Zcrf(q)W + (q)K-(-q) + ^(q)n + (-q) + n (q)n%-q) > 

T = 1 

\KK >= -^f(q)\K + (q)K-(-q)-K°(q)K°(-q) > (17) 

\irr ] >=Z<ff(q)W (q)r ] (-q)> 

where q is the momentum of the particles in the CM of the pair and q — \q\. 
We have used the convention that \ir + >= — |1, 1 > and \K + >= — ||, \ > isospin 
states. Note that for symmetry reasons there is no nir S-wave, T = 1 state. The 
function f(q) is normalized such that J2qf 2 (q) = 1 and we take it as infinitely peaked 
around a certain value q. Note also the apparent extra normalization factor l/y/2 
in the \tttt, T = > state which is a consequence of its symmetry (particles and 
antiparticles go into the same multiplet of isospin). By taking into account eqs. 
(17), (18) and the amplitudes in eqs. (3) to (15) we can write now: 
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V n = - < KK\£ 2 \KK >= -^(3s + Am], - £,p?) 

V 21 = - < 7rn\C 2 \KK >= -j^^s + 3m 2 K + 3m\ - § EiP-) (18) 
V 22 = - < ttttIAItttt >= (9s + ^ - 3 

T = 1 

V n = - < KK\C 2 \KK >= -^(3s - £ jft 2 + 4tt&) 

V 21 = - < n° V \£ 2 \KK >= ^0(6s -2E lP * + \m\ - \m\) (19) 

V 22 = - < 7i r]\C 2 \n°r] >= -^m\ 

where p\ = k , p 2 = p , = k' , = p' and the sum over momenta squared runs from 
1 to 4. For on shell amplitudes pf = mf. 

3 L = 0, T = 0, 1 strong amplitudes in the coupled 
channel approach 

It is easy to see that if one iterates the lowest order amplitudes calculated before, 
introducing loops, and one evaluates the finite contributions to the imaginary part of 
the amplitudes, the higher order contributions at energies around 1 GeV or below are 
even larger than those of the lowest order. Furthermore, we have the resonances fo ,ao 
in the L = channel around the KK threshold which should appear as singularities in 
the t matrix. It is then clear that the standard xPT expansion, keeping a few orders, 
neither will converge nor give this analytical structure. This difficulty is expectable 
since one is making an expansion in powers of p 2 /A^, and we go up to energies 
^/s = 1.2 GeV. In order to overcome these problems we assume that the lowest 
order Hamiltonian provides us with the potential that we iterate in the Lippmann 
Schwinger equation with two coupled channels, using relativistic meson propagators 
in the intermediate states. This idea follows identical assumptions made in refs. 



34] , |35H in the meson-baryon sector. Our channels are labelled 1 for the KK and 2 
for the 7T7T states in T = and 1 for KK, 2 for irr] in T = 1. The coupled channel 
equations then become 

t\i = Vn + VuGntxi + V\ 2 G 22 t 2 \ 

t 2 i — V 2 \ + V 2 \G\\t\x + V 22 G 22 t 2 i ^20) 

t 22 = V 22 + V 2 \Giiti 2 + ^22^*22^22 
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— mfj + ie (P — q) 2 — + le 

where P is the total fourmomentum of the meson-meson systems and q the fourmo- 
mentum of one intermediate meson. The terms VGT in eq. (20) actually mean 

VGT = J-0^V(k,p; q)G(P, q)t(q; k', p') (22) 

Note that as Ga has no angular dependence and is purely S-wave then ty is 
S-wave. 

The diagrammatic meaning of eqs. (20) is shown in fig. 1. 

The loop integral in eq. (22) is divergent. In xPT these divergences are cancelled 
by counterterms of chiral Lagrangians at higher order and some finite contribution 
remains from the loops and the counterterms. In our approach we take a cut off q ma x 
for the maximum value of the modulus of the momentum q. We vary this parameter 
in order to fix some property of the data (for example the position of the maximum 
of an experimental cross sections) and after this is done there is no freedom in our 
model. On grounds of chiral symmetry computations we should expect this cut off 



to be around 1 GeV 36 



By fixing the cut off one generates the counterterms in our scheme. It is also 
relevant to see that since our scheme produces the a , fo ,a resonances, our work 
makes connection with the work of ref. |37] where the finite parts of the counterterms 



are generated from the exchange of resonances. 

It is also clear that our non perturbative method does not generate all the terms 
which appear in the standard X-PT expansion. For instance loops in the t-channel 
are not accounted for. But the fact that the scheme proves so successful, as we shall 
see, indicates that one is summing the relevant infinite series needed to produce the 
structure of the t matrix. One can expect this because the singularities associated 
to these t-channel terms are far away from the physical region, contrary to what 
happens with the s-channel terms summed up in the Lippmann Schwinger equation, 
which lead to the pole structure of the amplitudes, essential to reproduce the physical 
scattering amplitudes. 

In principle one would have to solve the integrals in eqs. (20) by taking V and t 
off shell. However, this is not the case, as we show below, at least when dealing with 
S-wave, and we only need the on shell information. The argument goes as follows: 
As we can see from eqs. (18), (19), the on shell amplitudes are obtained by taking 
p 2 = m 2 and then we can write the off shell amplitudes as 



V = V on + (3j2(pt 



In order to illustrate the procedure let us simplify to one loop and one channel 
(the procedure is easily generalized to two channels). Hence we have 

V 2 = V 2 n + 2(3V on £(p? - m 2 ) + (3 2 $>, 2 - m 2 )(p 2 - m 2 ) (23) 

i ij 



7 



When performing the q° integration in the loop we have two poles, one for q° = 
Wi(q) and the other one for q° = P° + W2{q), where the indices 1 and 2 indicate the 
two mesons inside the loop. Let us take the contribution from the first pole (the 
procedure follows analogously for the second pole). From the second term in eq. (23) 
we get the contribution 

2(5V on r d 3 g (P»- Wi y- W l (3V on r 

(2vr) 3 J 2 Wl {q) {P* - w x ) 2 - w 2 2 (2tt) 3 J 19 1 ' 

which for large A compared to the masses goes as V on A 2 and has the same structure 
in the dynamical variables as the tree diagram. The third term in eq. (23) gives rise 
to the integral 

P 2 r d 3 q „,, p r s q 



(27r) 3 7 2w 1 (q) lK ' 21 (2tt) 3 J 2 Wl (q) 

coming from the first pole (q° = Wi(q)) and 



(5 2 r d 3 q 
J^Kfl 2w 2 (q) 1 

coming from the second pole (g° = P° + W2(q))- As we can see, the terms linear 
in P° cancel exactly, respecting the chiral structure which does not allow linear terms 
in P°. 

The term proportional to P 02 + m\ — m 2 in (25) leads again to a structure of 
the type [P 02 + (m\ — rn^)]K 2 and similarly happens with the quadratic contribution 
form the second pole which leads to [P 02 + (m 2 , — m^)]A 2 . These terms, together the 
V on A 2 which we obtained before, combine with the tree level contribution, giving rise 
to an amplitude with the same structure as the tree level one but with renormalized 
parameters / and masses. However, since we are taking physical values for /(/ = 
f n = 93 MeV) and the masses in the potential, these terms should be omitted. One 
can proceed like that to higher orders with the same conclusions. 

Since we are taking V and t on shells they factorize outside the q integral. Thus 
the term VGT of eq. (22), after the q° integration is performed by choosing the 
contour in the lower half of the complex plane, is given by 



VijGjjtjk — Vij(s)tjk(s)Gjj(s) 

(26) 

— JO (2^p"w 1 w 2 [P 02 -(^i+<-'2) 2 +j<:] 



where tUi = (q 2 + rri 2 ) 1/2 and P 02 = s and the subindex i stands for the two interme- 
diate mesons of the j channel. 

Thus the coupled channel Lippmann Schwinger equations get reduced to a set of 
algebraic equations: 

At = V (27) 
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where 











-I 




1 -I 






V h2 J 




V v 22 ) 



(28) 



1 — Viitu — V12G22 
^4 = I — V21G11 1 — V22G22 

-V21G11 1 - 

Introducing the notations 

A w = 1 — ^22^22 

A K = 1- VnGn 




A c = A K A W - V&GnGn (29) 
A = detA = A n A c 



and inverting the matrix A we obtain the following equations 

i 11 = ^-(AnVu + V&Gn) 
t 2 i = i : {V 21 G ll V ll + A k V 21 ) 

f — V22 1 y i2 G H 
r 22 — a r a a 



(30) 



where we have used the fact that V12 = V21 and ti2 = t2i by time reversal invariance. 

The physical amplitudes as a function of s (real variable) are given by the ex- 
pressions in eqs. (30). However, in order to explore the position of the poles of the 
scattering amplitudes one must take into account the analytical structure of these 
amplitudes in the different Riemann sheets. These sheets appear because of the cuts 
related to the opening of thresholds in Gjj(s) (see fig. 2). By denoting pi for the CM 
trimomentum of the K in the KK system and p 2 for the it in the 7nr, T = system 
or the 7r in the iri], T = 1 systems 

_ [s - (mij + m 2i f] l l 2 [s - (mii - m 2 Q 2 ] 1/2 

The sheets which we consider are 

Sheet I : Impi > 0, Imp 2 > 

Sheet II : lmp\ > 0, Imp 2 < 

Sheet III : Imp x < 0, Imp 2 < (32) 

Sheet IV: Imp 1 < 0, Imp 2 > 
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In order to make an analytical extrapolation to the II, III and IV Riemann sheets 
we have to cross the cuts Ga(P°). In order to do this we make use of the continuity 
property (for real P° > run + m<ii) 

Gf{P° + ie)=Gf{P G -ie) (33) 

where the index (b) indicates that we are in the second Riemann sheet of Ga while 
the index (a) indicates the first Riemann sheet. Then 



Gf (P° + ie) = Gf (P° - ie) = G^ (P° + ie) - 2UwGf (P° + ie) = 

(34) 

_ r {a),pO , - A , j \P 02 -(m li +m2 i ) 2 ] 1 / 2 \P 02 -(m li -m 2 i) 2 ] 1 ^ 
— ^ii V-r i* Lt ) t" 47r 2P 02 

We use the first and last terms of these equations to extrapolate G^ in the 
complex P° plane (the square root is taken such that Im yj > 0). Thus the sheet I 
is obtained with G$,G { £; the sheet II with G$l,G$ ; the sheet III with G$,G$ 
and the sheet IV with G§, Gf}. 



4 Results 



4.1 Phase shifts and Inelasticities: 

It is interesting to stress that we have only one free parameter at our disposal, which 
for chiral symmetry reasons should be of the the order of the 1 GeV. We have taken 
as our ultimate choice q max = 1.1 GeV which gives a good agreement with experiment 
in all channels for phase shifts and decay widths of the resonances. In order to obtain 



the phase shifts and inelasticities we use the two-channel S matrix [15 



S 



r/e 



2i<5i 



i(l-77 2 ) 1 /V(* 1+ *° Ve 



2iS 2 



(35) 



where 81,82 are the phase shifts for the 1 and 2 channels (KK,ttit in T = and 
KK, irt] in T = 1 respectively) and 77 is the inelasticity. The elements in the S matrix 
are related to our amplitudes via: 



£22 



2ipf( Sl1 ~ X ) 
•|f(522-l) 



(36) 



t 



12 



'21 



■^L(s 12 - 1) 



It is interesting to recall that the first two equations allow us to determine 81,82 
and 7] while the third equation allows us to determine 8% and 77. We obtain the same 
results with both methods, which is a check of consistency of the exact unitarity 
implemented in our scheme. In fig. 3 we show the results for the tot phase shifts 
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(Sow) in L = 0, T = 0. We show in the figure the results for two values of the cut 
off energy A = 1.1 and 1.2 GeV ( A = {q 2 max + m 2 K ) 1 / 2 ) 

We can see that the agreement with the experiment [1,38-45] is rather good up to 
values of \fs = 1.2 GeV . From there on discrepancies with experiment start to appear, 
but this should be expected because of the opening of the r\r\ channel, the increasing 
importance of the 4tt channel and the possible influence of higher resonances which 
cannot be considered as meson-meson states. We can see that the position of the /o 
resonance is well reproduced. On the other hand the results are moderately dependent 
on the cut off. At low energies, \fs < 500 MeV we can see that our results are rather 
independent of the cut off and agree with the data of ref. [fOJ. Another feature 
worth mentioning is a bump in our results around y/s = 500 — 800 MeV and in some 
experimental analyses, which in our case is originated, as we shall see, by the presence 
of the a pole. 

In fig. 4 we show our results for the phase shift of KK — > Tin, [8\ + 62) obtained 
through eqs.(35), compared to the experimental results of ref. [ iS| , [|7|] . The agreement 



with the data is also good. We have also calculated the results with A = 1.1 GeV 
and 1.2 GeV and although both are compatible with the data, we can see a better 
agreement with the data with smaller error for A = 1.2 GeV. The agreement with 
these two sets of data guarantees agreement with the KK phase shifts, which are 
obtained from these data by using the two channel unitarity of eq. (35). 

In fig. 5 we show the results for the inelasticity using A = 1.2 GeV. There is an 
obvious discrepancy between different data [45-48] . Nevertheless, we observe a good 
agreement of our results with the data of ref. [Q in consistency with the agreement 
found before for the phase shifts Sq KKtt7T coming from the same analysis pBj. It 



should be noted here that in ref. |L6| one obtains agreement with the upper set of 
data from ref. HH . 

For the T = 1 channel the data are scarce and different analyses on the same mag- 
nitude, which would show independence on assumptions made, are lacking. We can 
compare our results with data for irrj invariant mass distributions coming from the re- 
action K~p — ► S + (1385)7r~r7 [fE|. Similarly we can analyse K~p — > Y, + (1385)K~ K° . 
Following ref. [0 we write 

^ = C\t a \ 2 q a (37) 
dm a 

where m a is the invariant mass of the 7r~ 'rj(K~ K°) system, q a is the tt~(K~) momen- 
tum in the 7r~i](K~ K°) CM frame, t a is the scattering amplitude of the ir~r](K~K ) 
system and C a normalization constant. 

In fig. 6 we can see the n~rj invariant mass distribution and our results with two 
different normalization constants. In fig. 7 we show the same results for the K~K° 
invariant mass with the corresponding normalization constants. 

The peak in the 7i~r] mass distribution is due to the dominance of the ao resonance. 
We can see that we also produce the peak at the same place as the experiment 
and reproduce qualitatively the data. The agreement found with the K~K° mass 
distribution is also qualitatively good. 
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4.2 Poles and widths 



In order to determine the "true" mass and width of the a, fo and ao resonances 
one has to look for the position of the poles in the physical amplitudes in the complex 
P° plane. In the T = channel we find a pole in the II sheet of the complex plane 
which corresponds to the fo resonance at P° = (981.4 + i 22.4) MeV, corresponding 
to a mass Mf = 981.4 MeV and a width Tf = 44.8 MeV. This singularity results 
in a peak in the physical amplitudes tctc —>■ KK as one can see in fig. 8. 

We find another pole in the T = channel corresponding to the a. Actually we 
find two poles corresponding to zeros of A n and A c in eq. (29) which are very close. 
The zero of A n appears at P° = (468.5 + i 193.6) MeV while the zero of A c appears 
at P° = (469.5 + i 178.6) MeV. The fact that the two poles are so close and the 
associated width so large makes that in practice one only sees one bump in the cross 
section. In any case the second term of t 2 2 in eq. (30) is very small in the a region 
in the physical sheet (sheet I). It should be noted that the a pole in A c also affects 
the in, £21 amplitudes, but this has scarce practical consequences since this appears 
in the unphysical region below threshold and the physical region near threshold is 
dominated by the f pole. 

The a that we find is rather broad with IV ~ 400 MeV. We should also note 
that since it comes essentially from the first term of £22 in eq. (30), it is practically 
tied exclusively to the pion sector, with practically no influence from the KK sector. 
Hence, if we eliminate the KK channel from our coupled equations the a pole survives. 

In the case of the fo we have also observed that it is dominated by the KK channel, 
in the sense that if we eliminate the tttt channel (setting V12 = 0) the fo pole appears 
also, now at P9 = (973.4 + iO) MeV. The zero width is obvious in this case since 
there is no decay into irn, nor possible KK decay because one is below threshold. In 
the case when the tttt decay channel is open, the KK channel also opens because the 
finite fo width gives phase space for the decay into KK (we will come back to this). 

In the T = 1 channel we find a pole in the II sheet corresponding to the ao at a 
position P° Q = (973.2 + ^85.0) MeV . This corresponds to a width of T ao = 170 MeV. 
This pole is much tied to the coupled channel. Indeed, if we eliminate the f]ir channel 
(setting Vu = 0) the pole disappears. 

In fig. 8 and 9 we show the results for the |ti2| 2 around the f and a regions 
respectively. We see there two peaks corresponding to the resonance excitation around 
y/s = 987 MeV for the fo and 982 MeV for the ao- The naive widths that one 
observes there are around 20 MeV and 60 MeV respectively. The position of the 
peaks and the "widths" do not correspond exactly to the numbers of the pole position. 
This is a well know fact of resonances. However, in the present case there is one extra 
feature, which is mostly responsible for the shrinking of the width and change of 
position of the maximum. This was discussed in ref. [OT for the specific case of 
the ao, but the arguments hold also for the fo. This feature is a cusp effect, due to 
the strong coupling of the resonances to the KK channel which appears around the 
pole position. This makes the KK decay width to increase very fast as the energy 
increases, reducing the amplitudes accordingly and producing a shift of the peak 
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position. Large widths for the ao around 200 MeV have also been claimed in former 
studies IlIB] . 



5 Strong partial decay widths 

The particular structure of the amplitudes in, £21,^22 of eq.(30) makes that the tu 
and t\2 amplitudes in T = have a peak around the fo pole, while the tttt amplitude 
combines the fo pole with a background coming from the terms that give rise to the 
a pole and one has actually a minimum around that region. In the case of T = 1 all 
these amplitudes peak around the a pole. In view of this, we take the tu and £ 12 
amplitudes in both cases and parametrize them around the resonance region in terms 
of a Breit Wigner expression with energy dependent widths. 

y. _ 1 g|l(f) 

tn - 2M R ^_ MR+i EiM 

f _ 1 gfii(*Wa) ( 38 ) 
121 - 2M R ^_ MR+i EiM 

We shall not need the explicit values of gm(s) and T R (s) since we can compute 
our final expressions in terms of lmtj\ as shown below. 

In order to evaluate the partial decays widths we take into account the finite width 
of the resonance and substitute in the standard formula for the decay width 

2w5(M a -ui- cu 2 ) 21m — ^- (39) 

Mr- UJ X - UJ 2 - l^t 

in view of which the partial decay widths are give by 

r - = is? L +m dW w b»w\ \ MR _ w / +m («) 

However in practical terms we write eq. (40) in terms of the amplitudes of eq. 
(36) directly and we find 

IV = C 1+m2 dW 4,AM R Imt xx 

(41) 

r«,2 = -jh Q +m > dW &4M R ®£>- 

We integrate in W around two widths up and down the peak (if allowed by the 
lower limit) where we find that Tr^+Tr^ = T tat , where Y tot is the width corresponding 
to the pole of the amplitude evaluated before. 

In table I we give the results of the partial widths and the mass of the /q. Both 
the pole position and the peak of the t 2 ± amplitude are in agreement with the results 
of the particle data group ]5T| . As with respect to the width, Y tot is also compatible 



with most of the present data as can be seen in ref. | 51fl . The ratio T^/iT^ + T K ^) 
is also in good agreement with the different measurements, as can be seen in table I. 
Since KK branching ratio is obtained by subtraction from unity of the itit branching 
ratio, our results are also compatible with those quoted in the particle data book. 
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In table II we give the results for the Oq. The mass quoted in the PDB [51j is in 
agreement with our results for the peak position of the £21 amplitude. Although the 
position of the pole of the amplitude appears at a different energy. This difference 
between the pole position and the peak of physical amplitudes is well known in general 
g|. The width of the resonance is Ytot — 170 MeV, apparently larger than many 



of the experimental numbers in []51j] . However, we already quoted how the cusp 
effect can convert this large width into experimental narrower amplitudes |5(| and 
this was indeed the case in our results as shown in fig. 9. We also quoted our 
qualitative agreement with other theoretical results JTTJ which have r ao = 202 MeV. 
Another interesting piece of data is the fraction of KK decay into tjtt. This fraction is 
practically unity in our approach and compatible with the most recent measurement 
of ref . . We should note that the fact that the fraction of KK to t/tt decay is large 
is due to the fact that the width is large, since for narrow resonances there would 
be little phase space for KK decay. 

As we see, we can claim a global good agreement with experiment for masses and 
partial strong decay widths. 



6 Decay of the resonances into 77 

In order to determine the decay width of the cto and fo resonances into the 77 channel 
we proceed in an analogous way as done for the strong decays. In the first place we 
obtain the KK — > 77 amplitude using the coupled channel approach discussed above, 
and second we obtain the effective coupling of the resonances to the 77 channel in 
order to obtain the partial decay width via eq. (41). The Lippmann Schwinger 
equation for the amplitude KK — > 77 would now be written as 



tyl — Vyl + V^lGntu + ^2^22^21 



(42) 



where now V y i, V l2 are the potentials obtained from the chiral Lagrangian at lowest 
order with external electromagnetic currents 



V- 



■yyK+K- 



-2*e 2 {eh 



k + -k 2 
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(43) 



7777 ^ 



— lie 2 le^e^ — p +' eip -' 62 _ p+-^p--^i \ 



64| we have k\ , k 2 the photon 
the momenta of the 



In eq. (43), following the nomenclature of ref. 
momenta, ei,e 2 their polarization vectors, and k + , k^,p + ,p 
K + , K~ , 7r + , n~ respectively. We also take the same gauge as in ref. |63| in order to 
perform the calculations, ei • £7 = ei • k 2 = e 2 • A7 = e 2 • k 2 = 0. From eqs. (43) and 
the isospin states of eqs. (16), (17) we obtain the S-wave V 7 i potentials. We factorize 
the on shell t matrices and Vyi outside the integrals in VyiGuTu in eq. (41), as we 
did in the strong amplitudes, and we evaluate the integrals of Gu cutting the integral 
over momentum at the same value of q ma x- Since the scattering amplitudes tn are 
evaluated before, the evaluation of t 7 i of eq. (42) is straightforward. Since both t±i 
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and t 2 i become singular at the f , a poles, the t 7 i amplitude also has a singularity 
there in the T = 0, T = 1 channels respectively. 

It is interesting to recall than the standard xPT expansion acounts for some terms 
not accounted for here pi. Once again, if these terms were added to our amplitude 



they could not change the resonant part of the t 7 i amplitude of eq.(42) given by tu 
and £21 and which is observed experimentally for 77 — > 7r + 7r~, 7r°7r°, 7c°r] 
By following similar steps as for the strong decay we can write 



1 1 f°° q (Imt 



2 167T 2 J-mi+m 2 VT 2 Imtn 

where mi, 777.2 refer now to 2m. K , [m^ + m v ) for T = 0, (T = 1) since below this value 
there is no imaginary part of JmT l7 . The factor 1/2 in eq. (44) is due to identity 
of the two final photons. As we did before, we integrate over W in an interval of 2T 
above Mr and 2T below, if allowed by phase space. 

The results for the widths are shown in tables I and II for the /o an d clq respectively. 
We see a very good agreement with the PDB averages which are compatible with many 
of the different results of the analyses quoted in those tables. In any case the situation 
is still problematic, as it is reflected in table I in the dispersion of the quoted results 
of the analyses. In fact in this sector,T = and L — 0, it is not understood why a 
broad scalar structure above lGeV is observed for 77 — > 7r + 7r~, following ref. [4], 
and it is not seen in 77 — > 7r°7r°. On the other hand for the T = 1, L = channel the 
problems come from the way one isolates the background around the well observed 
do (980) and 02 (1320) resonances in the 77 — > n t] reaction, with its influence on the 
coupling to two photons of the former resonances. All these problems are properly 
reviewed in [f26fl. 
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7 Conclusions 

We have used a non perturbative approach to deal with the meson meson interac- 
tion in the scalar sector at energies below \/s ~ 1.2 GeV, exploiting chiral symme- 
try and unitarity in the coupled channels. The Lippmann Schwinger equation with 
coupled channels, using relativistic meson propagators, and the lowest order chiral 
Lagrangians, providing the meson-meson potential, are the two ingredients of the 
theory. One cut off in momentum, expected to be of the order of 1 GeV, is also used 
in order to cut the loop integrals and it is fine tuned in order to obtain essentially 
the position of the poles. This is the only parameter of the theory. The best fit to 
all the data is found with q max ~ 1.1 GeV, or equivalently A = 1.2 GeV, both for 
the T = and T = 1 channels. We find poles for the fo and ao resonances, for T = 
and T = 1 respectively, and also for the a in the T = channel which appears as a 
broad resonance (r = 400 MeV) at about 500 MeV of energy. 

We compute the T = phase shifts and inelasticities for mi — » tttt, tttt — > KK and 
KK —>■ KK and they are in good agreement with experiment. In the case T = 1 we 
look for mass distributions of ti Q 7] and KK which are also in agreement with poorer 
experimental data. 
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We also compute the partial decay widths of the irir, KK and Tc r], KK respectively 
and the results are in good agreement with experiment. We have also studied the 
partial decay mode of the resonances into 77 . We show that the KK — > 77 amplitude 
(as well as other ones) contains the f and ao poles in the T = and T = 1 channels 
respectively and we have evaluated the decay width of the ao and fo resonances into 
the 77 channel. The results obtained are in good agreement with experiment. 

The scheme used has proved to be very successful. The amount of data that one 
can reproduce using one only parameter is really impressive. We should recall that 
fits to some of the data studied here required a large number of parameters. On 
the other hand it is remarkable to see that this parameter is actually of the order of 
magnitude of what we could expect from former studies of chiral perturbation theory. 

The teaching of the present results is that the constraints of chiral symmetry at 
low energies and the implementation of unitarity with the coupled channels, keeping 
also the real part of the loops, is the basic information which is contained in the 
meson-meson interaction in the scalar sector below y/s = 1.2 GeV. 

On the other hand it is quite clear that the present scheme has allowed us to make 
progress within the lines of chiral dynamics, introducing a non perturbative method 
which allow us to go to higher energies than in xPT and at a lower price in the 
number of parameters. It becomes appealing to exploit these ideas in other processes 
related to the scalar sector or try to extend it to other sectors in order to test the 
potentiality of the method and its limits. 
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Table I: f Mass and partial widths (A = 1.2 GeV) 



fo 


our results 


experiment 


Mass 






Pole position 

[MeV] 


981.4 


980 ± 10 [51] 


Mass 






Peak of the t^i 
amplitude 

[MeV \ 


987 


980 ± 10 [51] 


Ttot [MeV] 


44.8 


40 - 100 [51] 






0.67 ± 0.009 [52] 


^ TV TV 


0.75 


o.8i ±8:gg [53] 

0.78 ± 0.003 [54] 
(av.) 0.781 ± 0.024 [5lJ 






0.63 ± 0.14 [4] 


r 




0.42 ±0.06 ±0.18 [55] 


i 77 

(KeV) 


0.54 


0.29 ±0.07 ±0.12 [56] 




0.31 ±0.14 ±0.09 [57] 
(av) 0.56 ± 0.11 [51] 
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Table 2: a Mass and partial width (A = 1.2 GeV) 



a 


our results 


experiment 


Mass 






Pole position 

[MeK J 


937.2 


983.5 ± 0.9 [51] 


Mass 






reak ol thei2i 
amplitude 

[MeV] 


982 


983 ± 0.9 [51] 


rtot [MeV] 


170 W 


50 - 100 [51] 




1.07 


1.16 ± 0.18 [58] 
0.7 ±0.3 [59] 






0.25 ± 0.008 [60] 




0.22 


0.28 ±0.04 ±0.01 [55] 
0.19 ± 0.07 ±g ? 7 [61] 
av.0.24±°;™"[51] 



(*) See discussion in the text about the narrower physical amplitudes because of the 
cusp effect. 
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FIGURES 

Fig.l. Diagrammatic series for eqs. (20) for with T = 

Fig. 2. Sheets in the complex energy plane. The unitarity cuts are indicated by 
thick lines. 

Fig. 3. phase shifts for kk — > kk in T = and L = compared with different 
analyses indicated between brackets. 

Fig. 4. ^qkRittt P nase shifts for KK — > kk in T=0 and L=0 compared with 
different analyses indicated between brackets. 

Fig. 5. in T = and L = compared with different analyses indicated 

between brackets. 

Fig. 6. k~t] mass distributions. The data are from the reaction K~p — > T, + (1385)k~ i] 
[49]. 

Fig. 7. K~K° mass distributions. The data are from the reaction K~p — > 
£+(1385)^-^° [49]. 

Fig.8. |t; 2 | 2 , T = 0, around the mass of the / (980). 

Fig.9. \t\ 2 \ 2 , T = 1, around the mass of the ao(980). 
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